typeset using JPSJ.sty <ver.l.Ob> 



o 
o 
o 

(N 
D 
00 



C/3 



Numerical Study on Aging Dynamics in Ising Spin-Glass Models: 
Temperature- Change Protocols 
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By means of Monte Carlo simulations on the three-dimensional Ising spin-glass model, we have 
studied aging phenomena with various temperature(r)-change protocols. Particularly, a T-shift 
protocol, in which a system is first quenched to and aged for a period twi at a temperature Ti, 
and subsequently aged at a new temperature T2 is closely investigated. Most importantly, the 
mean size of domains, which is extracted from the replica-overlap function, was found to grow 
monotonically without any appreciable decrease by the T-change. We also found the relaxation 
of energy density and spin auto-correlation function during the T-shift process can be explained 
within the following picture: the dynamics finally crossovers to an isothermal aging of T2 at 



around £2 — after the T-change, where t)^^' is an effective waiting time related with iwi, 
but in the transient regime £2 ^ ^wf "* > adjustment of the population of thermally active droplets 
from that at Ti to T2 takes place very slowly. Implications of the results on aging phenomena in 
other T-change protocols such as T-cycling and continuous T-change with an intermittent stop 
observed by simulations as well as experiments are discussed. 
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§1. Introduction 

Aging phenomena in spin glasses have been extensively 
studied in recent years. Experimentally, many elabo- 
rated protocols have beeii.^adopted to reveal various as- 
pects of the phenomena.EfH'H^ Theoretically, on the other 
hand, many ideas and models on spin glasses and re- 
lated systems have been proposed and examined to get 
a proper_imderstanding of aging dynamics in such glassy 
systems.QB On the side of numerical study, a number of 
simulations, which could account for qualitatiyje features 
of the aging effects, have been performed.ErQ'H-' In spite 
of the enormous efforts, however, a unified picture of ag- 
ing dynamics in spin glasses has not been established 
yet. The difhculty is directly connected to nature of the 
low-temperature spin-glass (SG) phase which has been a 
controversial issue since more than a decade ago. 

Among simulational results on the short-ranged Ising 
SG models, of particular importance is the growth of a 
mean domain size which is observed through the replica- 
overlap- Eunclion, G{r,t), in an isothermal aging pro- 
cess. Bllj^llil't3 In the process a system is quenched in- 
stantaneously from temperature above its SG transition 
temperature to T below Tc, and is isothermally aged 
at T. The function G{r, t) is defined as 



G(r,i) 



1 ^ 



:(/3) 



(1.1) 
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where a and P are indices of two replicas which have dif- 
ferent random spin configurations at i = (an instant of 
the quench) , and are independently updated. Its correla- 
tion length Rt (t) is found to grow almost in a^pojver-law 
in the case of 3 dimensional (3D) modclsJ13'l 



(1.2) 



where Lq and tq are certain characteristic units of length 
and time, respectively, and the exponent 1/ z{T) depends 
on T. In our previous workJlS* hereafter referred to as I, 
we obtained 



l/z{T) ~ 0.17T, 



(1.3) 



at T < 0.7 J for the Gaussian model with zero mean and 
variance of interactions J. It is also fitted to a logarith- 
mic functional form written as i?Wt)-pi,i?o + b{\nty/^ 
with Ro,b and -0 being constantsB't2lt3 This ambigu- 
ity is attributed to the insufficiency of the time window 
of simulations which can be performed by the presently 
available computers. Although there remains such an 
ambiguity in the explicit growth form of Rxit), we can 
naturally regard Rt (t) as a characteristic length scale of 
the SG ordering in aging dynamics, i.e., a mean distance 
of domain walls separating different pure states. 

It is natural to ask to what extent the mean domain 
size Rrit) is related to aging behaviour obser^d in yaxj- 
ious quantities simulated. The droplet theory,liSE3'llj' \3 
particularly the one proposed by Fisher and IIuse,L3' 
contains scaling ansatz which relate various quantities 
to the characteristic length scales such as Rxit). We 
have been numerically studying aging phenomena in 
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SG model to clarify the validity of the 
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Fig. 1. Schedules of T-change in some experimental protocols of 
aging: 1) isothermal, 2) T-shift, 3) T-cycle, and 4) continuous 
T-change with an intermittent stop. 



Experimentally, various interesting phenomena have 
been also observed in aging processes in addition to the 
isothermal one. They are T-shift, T-cycling, and con- 
tinuous T-change with an intermittent stop(s) protocols 
whose schedules of T-change are schematically shown in 
Fig. 1. The ac susceptibility has been continuously mea- 
sured through these T-changes. In the case of a negative 
T-shift process with AT = T2 - Ti < in Fig. 1 and 
with a magnitude |AT| larger than a certain value, the 
ac susceptibility is observed to increase discpntixiiiDusly 
just after t = t^i and then to relax againJlj'EHlLlI' Its 
behaviour looks as if the system restarts aging by the 
T-shift at t = iwi forgetting its thermal history before 
iwi- It is called the rejuvenation, or chaos effect. Fur- 
thermore when the temperature is turned back to Ti at 
t = twi+iw2 (T-cycling protocol), the value of the ac sus- 
ceptibility returns to that on the isothermal aging curve 
of Ti at almost t — i^i, and then relaxes along the curve. 
This behaviour indicates that the system in fact remem- 
bers its history before and is called the memory 
effect. Similar coexistence of the paradoxical phenom- 
ena, rejuvenation and memory effect, has been recently 
observed also in-a-continuous T-change with an inter- 
mittent stop(s).c3'E3 These experimental findings have 
stimulated theoretical interest .EP 

The main purpose of the present paper is to investi- 
gate aging phenomena in various T-change protomls by 
simulations on the 3D Gaussian Ising SG model.cS We 
concentrated mostly on the T-shift process, which is the 
most fundamental among the various T-change proto- 
cols, since the latter are considered as certain combina- 
tions of the former. We investigated the time evolution 
of the mean domain size and the relaxation of the energy 
density and the spin auto-correlation function during the 
T-shift process. 

We have found, within the time window of our simu- 
lations, the following characteristic features of the aging 
dynamics, i) The mean domain size R[T2,Ti] ^wi) in the 
T-shift process (see Fig. 1), with either a negative or pos- 
itive AT, does continue to grow. Here and hereafter t is 
a time measured from the first quench from above Tc to 
Ti, t2 — t — iwi a period that the system ages at T — T2 
after the T-shift a,t t — t^i- ii) The time evolution of 



the energy density and the spin auto-correlation func- 
tion is related with the mean domain size -RjjVj j-^j (t; t^i) 
at sufficiently long time scales after the T-shift by the 
scaling relations which are the same as isothermal aging, 
iii) However, there is a transient time regime after the 
T-shift where there is extra contribution to the relax- 
ations. Interestingly, the extra contribution also shows 
very slow, non-exponential relaxation. We interpret the 
latter as not djj£_.to the chaos effect predicted by the 
droplet theorytSllIP but as due to adjustment of the pop- 
ulation of droplet excitations within each domain. 

The memory effect observed experimentally is easily 
interpreted from the above-mentioned scenario, namely, 
it is attributed to the persistence of domains in a T-shift 
process. It is not yet conclusive, however, whether na- 
ture of the transient regime in the scenario is common to 
phenomena observed experimentally just after a T-shift. 
In this respect, the recent work on aging in the ferro- 
magnjctic fine particles system (FFFS) done by Mamiya 
et alES' is of quite interest, since the time window of their 
observation in_ujiit of the microscopic time ('^ 10~^sec 
for this systemcj) ) is rather closer to that of our simula- 
tions. Indeed, they have observed T-shift aging phenom- 
ena which can be well interpreted by our scenario. It 
seems, however, that the scenario is hard to explain the 
rejuvenation, or chaos-efiecl. pbserved experimentally in 
ordinary spin glassesJ13'EilE3'E§' whose microscopic time 
is of the order of 10"^'^ sec. 

The present paper is organized as follows. In the next 
section we briefly review the scaling properties of isother- 
mal aging derived from the droplet theory, add further 
comments on them, and present the scenario on the T- 
shift process. In §3 we present the results of our simula- 
tions on various aspects of the T-shift, as well as the T- 
cycling and continuous T-change processes. In the final 
section implications of the results with the experimental 
observations are discussed. 

§2. Phenomenological Picture 

2. 1 Scaling properties of isothermal aging 

Here we briefly review the scaling propertiasuaf isother- 
mal aging derived by the droplet thcoryllJ't3 due to 
Fisher and Huse (FH). In an isothermal aging process 
a system is quenched from above to temperature T 
below Tc. At waiting time t^ after the quench, there 
are domain walls separating different pure states of the 
SG phase at a typical distance i?T(^w) from each other. 
The distance Rxitw) can be measured in simulations for 
instance through the replica overlap function defined in 
eq.(y). 

Except for places where the domain walls run, the bulk 
of the system is essentially in equilibrium. An impor- 
tant subtle feature inside the domains of a spin glass is 
that there can be large scale thermal fluctuations due to 
droplet excitations. In ideal equilibrium, a droplet exci- 
tation can be defined as a global flip from a ground state 
of a droplet (cluster) oLspins within a distance L/2 from 
a certain given site iS3 The latter can be considered 
as a simple two-state system with a free-energy excita- 
tion gap FL(i) and a barrier energy BL(i) of a thermal 
activation process. 
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The typical value F^^ of the gap FL^i) scales as, 



(2.1) 



where T is the stiffness constant of domain walls on the 
boundary of the droplet. The gap is however broadly 
distributed and its distribution function pl{F) is con- 
sidered to follow the scaling form, 



Pl{F) 



1 



ptyp ' 
^ L 



(2.2) 



A very important property of the scaling function p{x) 
is that it has finite intensity at x = 0, p{0) > 0, which 
allows large scale droplet excitations even at very low 
temperatures. The probability to have a thermally active 
droplet at a very low temperature T is given by 



p{L\T) 



rP(0). 



(2.3) 



T(i/Lo)^' 

The latter feature yields equilibrium properties which 
make spixL-glasses distinctly different from simple ferro- 
magnets.E3^ 

The relaxation time of a droplet of size L centered at 
site i is given by the Arrhcnius law 



tl(i) = To exp 



fcRT 



(2.4) 



where B^^i) is the energy barrier of the droplet. The 
typical value Bf^''^ of the energy barrier -8^(1) scales as, 



i3fP~A(i/Lo)^ 



(2.5) 



where A is a characteristic unit for energy barriers as- 
sociated with thermal activation of droplet excitations, 
and the exponent ip satisfies 9 < ip < d — 1. 

In the droplet theory, it is assumed that the growth 
of domains is also governed by droplet excitati ons . An 
obvious but very important consequence of eqs.(2.4) and 
(2.5) is that dynamical processes, including both the do- 
main growth and droplet excitations, at different length 
scales are extremely widely separated.|j^ Hence within a 
time scale of small scale processes, large scale processes 
look as if they are almost frozen. 

In isothermal aging, for example, the relaxation of the 
excessive energy per spin Serit) with respect to the equi- 
librium value is expected to follow the scaling form as, 



Serit) f{RT{t)/Lof/{RT{t)/Lo)'', 



(2.6) 



where T is a characteristic energy scale, and d is the 
dimension of space which is 3 here. This scaling form 
was confirmed, in our previous work I, with the energy 
exponent 9 = 0.20 ± 0.03 which agrees with the result of 
the defect energy analysis at T = 0E3 

The droplet theory also predicts important scaling 
properties of the spin auto-correlation function. 



(2.7) 



in isothermal aging. Here the over-line denotes the av- 
erages over sites and different realizations of interactions 
(samples), and the bracket denotes the average over ther- 
mal noises. In the so-called quasi-equilibrium regime de- 
fined by r <C iw, the typical size Lt{t) of droplet exci- 



tations which take place in the time scale of t is much 
smaller than the typical separation Rxit^) of domain 
walls present after the waiting time For this situ- 
ation the droplet theory introduces a phenomenological 
concept called effective stiffness which characterizes a re- 
duction of the excitation gap of small scale droplets due 
to the presence of the domain walls. The latter yields 
tw-dependence of CxiT^t^) which is written in terms of 
Rritw) and Lt{t) as the following: 



l-CT(T;iw) = ACT(r;tw) +aT(T), 



with 



(2.8) 



(2.9) 



and 



ACTir;t^) ^ ml ^ ^° 



( Lt{t) 



'T \Lt{t)J \RT{t^), 

(2.10) 

Then the zero field cooled susceptibility Xzfc:t('''; tw) 
which is often observed by experiments should also be 
understood within the same ansatz since the fluctuation- 
dissipation theorem (FDT) 

XzFC;T(r;tw) - - CT{T;t^)] = XT{r;t^). (2.11) 

holds in the quasi-equilibrium regime. The FDT 
has bfien-|tested and confirmed by the careful experi- 
menti£*]lH) even in the presence of weak violation of the 
time translational invariance. 

2.2 Further comments on isothermal aging 

In our previous work,E3^ hereafter referred to as II, 
wc confirmed the scaling ansatz on CT(T;tw) described 
in the previous subsection. Here, we first make the com- 
ment that Xzfc;t("'" ; ^w) experimentally observed and an- 
alyzed through eq.( 2.11 ) follow the scaling ansatz. 

In Fig. 2, for example, we plot Xzfc;t('''; iw) of a 
Ag:Mn 2.6at% spin glass against t^/r. It is evaluated 
from totrm;t('''; iw) measured by the Saclay groupliP as 

XZFC;T(T;t„) = [tofC;T -TOTRM;T(r;tw)]//i with mFC;T 

being the field-cooled magnetization under a field of 
strength h. We note that we plot here only a part of 
their data with ^ t, i.e., those in the quasi-equilibrium 
regime. The data for different r's lie on a universal curve 
when they are vertically shifted by suitable amounts de- 
pending on T as shown in the inset of the figure. The 
feature is quite similar to 1 — CT(T;tw) analyzed in II. 
[Note that the factor {Lo/Lt(t))^ is practically con- 
stant in eq.( ^.lO ) since 9 is much smaller than d — 9, as 
was the case of 1 — Ct(t; t„) in II.] This indicates that 
Xzfc:t(t; t„), via eq.(|2.11), follows the scaling forms 



eqs.(^), (|2^ and (|2.10|) as well. The values of the 
exponents (d — 9)/z{T) and 9/z{T) turn out to roughly 
agree with the estimates obtained by making use of our 
previous results in II, 9 ~ 0.20 and l/*ff) ^ 0.16(T/rc), 
i.e., eq.(|L|) with Tc = 0.95 ± 0.04,B However, such 
semi-quantitative agreement is lost for the data at tem- 
peratures closer to T^. 

An intriguing point is that the above result sug- 
gests that Rxit) in the time scale of experiments is 
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Fig. 2. XZFC;T{T;tw) of a AgMn spin glassfl at T/Tc = 0.87 
for fixed r's plotted against t-w/r. Only the data in the 
quasi-equilibrium regime are plotted, namely, they are for 
T = 1, 2, 4, 8, 256sec from bottom to top with tw = 
300, 1000, 3000, 10000, 30000sec. In the inset the scaling be- 
haviour of Axzfc tC''"; iw) corresponding to eq.(2.1C) is demon- 
strated. 



also compatible with the estimate from our simulations. 
lUj-this context it is worth pointing out that Joh et 
alE3 have recently extracted Rrit) of Cu:Mn 6at% and 
CdCr1.7Ino.3S4 spin glasses from the thermoremanent 
magnetization (TRM). The growth law they have ob- 
tain ed t urns o ut to agree with the simulational result of 
eqs.(p^) and (L3) even quantitatively. The estimated 
Rrit) with t of the order of laboratory time is at most 
only several tens of lattice distances. However, the ex- 
perimental studies on the isothermal aging dynamics so 
far reported have not converged to a unified picture yet 
as we discussed in II. 

Another comment is related to the characteristic 



length and time scales, Lq and tq in eq.(1.2), which we 
did Jiipt fix explicitly in II. According to the droplet the- 
ory they are some microscopic scales at lower temper- 
atures, but near below Tc they are expected to crossover 
to the SG coherence length ^_ and the corresponding 
critical relaxation time, respectively. Indeed, Hukushima 
et alSy have found recently such a crossover behaviour 
in aging dynamics in the 4D ± J Ising SG model. A sim- 
ilar analysis on the present 3D model certainly has to be 
done, which may clarify the disagreement between the 
experimental and simulational results close to Tc men- 
tioned just above. Furthermore, as suggested in Ref. 34, 
it may also provide a useful key to resolve the apparent 
ambiguity of the growth law mentioned previously. 

2.3 T -shift process within overlap-length 

Now let us consider the T-shift process. After a system 
is first quenched from above Tc to a temperature Ti be- 
low Tc, domains of different pure states grow up. Their 
mean size reaches to i?Ti (^wi) after a given waiting time 
t = tvri- After the temperature is changed at t = t^i to 
a new temperature T2, we expect the domains continue 
to grow but with a rate specific to the new temperature 



T2 (see eq.( p.l[ ) below). It can be said that, at the in- 
stant of the T-shift, the system is aged as if it is aged 
for the effective waiting time by isothermal aging 



at T2, where t^f ^ is determined by 

Rxiitwl) — ^^Tsl^wf'') 



(2.12) 



with Rrit) being the growth law of the mean domain 
size in isothermal aging at T. 

Another important process which must take place after 
the T-shift is adjustment of the population of thermally 
active droplets within the domain. Since the equilibrium 
population is proportional to T as given in eq. ( ^.3[ ) , there 
must be a change of the population of an amount pro- 
portional to AT = T2 — Ti at each length scale L. It 
is emphasized that this process itself is intrinsically very 
slow since adjustment of the droplets takes place only 
by activation processes described by eqs.(^.4|) and (pTsj). 
Thus we naturally expect that there is a new charac- 
teristic length T[y2,Ti] (^2, ^wi) which slowly grows after 
the T-shift such that the population of droplets smaller 
than Tjj-^ (^2, ^wi) is equilibrated to the new tempera- 
ture T2. However, the population of droplets larger than 
L[T2,Ti]{t2,twi) is still equilibrated with respect to the 
old temperature Ti. We call the characteristic length 
-^[T2,Ti] (^2, iwi) a size of quasi- domains. 

The growth law of Tjj-^.Ti] (^2, ^wi) is expected to be 
the same as the domain. When the system ages at T2 for 



a period about t^°P defined by eq.( ^.12 ), Tj^^.Ti] (*2, ^wi) 



catches up R[T2.Ti]{titwi) and the quasi-domains merge 
into domains in isothermal aging at T2. Thus at large 
time scales t2 ^ i^f , aging becomes essentially the same 
as isothermal aging at T2. The time range ^2 ^ ^wf ^ on 
the other hand, is regarded as a transient regime where 
adjustment of the population of droplets are taking place. 

In the above argument we implicitly assumed that 
equilibrium states at different temperatures are the same 
except for the change of the population of droplets 
eq.(2.E). It should be remarked that the assumption 
makes sense only for lengtli_scales smaller than the so- 
called overlap length /AT-t^'Ea' It is predicted that equi- 
librium states at different nearby temperatures Ti and 
T2 differ from each other at length scales larger than 
this length Ia'p-i This aspect is called chaotic nature of 
the SG phascES) The overlap length I^t is expected to 
diverge when AT = T2 - Ti as 

2/(d^-2e) 



Iat ^ Lq 



T3/2 



Ti/2|AT| 



(2.13) 



where, 9 is the exponent in eq.(2.1) and ds the fractal 
dimension of surface of the droplets. 

§3. Results of Simulations 

3.1 Model and method 

We have carried out MC simulations on various ag- 
ing processes in the same 3D Ising SG model as in 
our previous works, i.e., the one with Gaussian nearest- 
neighbor interactions with mean zero and variance J = 1 . 
The critical— temperature has been obtained as Tc = 
0.95 ± 0.04.23' The heat-bath MC method we use here 
is also the same as in I and II. The data we will discuss 
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below are obtained mostly at T = 0.6 ^ 0.8 in systems 
with linear size Ls = 24 averaged over a few hundreds 
samples with one MC run for each sample. In I, it was 
confirmed that finite size effects do not appear for these 
values of T and Lg within our time window ( < 2 x 10^ 
MCS). 

3.2 SG coherence length — domain growth 

Let us begin with time evolution of R[T2,Ti]{'t':twi)i a 
mean domain size in a negative T-shift protocol with 
AT = T2-Ti < (see Fig. 1). The simulated results 
with Ti = 0.8, Ta = 0.6 and t^i = 1000 MCS are pre- 
sented by solid circles in Fig. 3(a). They are extracted 
from the replica overlap function cq.(y_T|) which is ob- 
tained by averaging over several thousands samples with 
L, - 32^ 

Clearly, i?^^^ 7-^] (i; t^i) does grow continuously with- 
out any noticeable decrease by the T-shift. This im- 
plies that, just after the T-shift, the system looks as 
if it has aged isothermally at T = T2 for a period of 



the effective waiting time t^°f^ defined in eq. ( pl^ ). In- 



deed, when R[T2.Ti]{t',tvji) at t > <wi is plotted against 
i = t^i^ + 12 (open triangles in the figure) , it almost lies 
on the isothermal aging curve (t), i.e., the two curves 
are related as 



R[T2,Ti]it;t^ 



^T.(i2 + 4f^), at t2 > 0. (3.1) 



Here, a value of i^f (~ 6230) has been chosen so as to 
obtain a good collapse of the two curves at time range 



<2 ^ ^wf The chos en v alue t urns out to satisfy eq.( 2.12 ) 



combined with eqs.(L2) and (1.3). 

As shown in Fig. 3(b), -RfT^.Til(^; ^wi) in a positive 
T-shift protocol is also described by eq.( t3.l[ ), but with 
4°f ^ (- 1400) < iwi (= 10^) in this case. The contin- 
uous growth of domains in the T-shift process, either a 
negative or positive one, is one of the most important 
results observed in the present simulation. 

3.3 Evolution of energy — Transient and isothermal 
regimes 

In Fig. 4(a) we present the time evolution of energy per 
spin, e[T2,Ti]itjtwi), in the same negative T-shift process 
as shown in Fig. 3(a). The original data drawn by the 
solid circles quickly deviate from the isothermal aging 
curve, eTi{t), just after the T-shift at t = i^i, and then 
slowly crossover to (t) from below at longer times. 
As drawn by the open triangles in the figure, the part of 
e[T2,Ti]i't': twi) after the crossover can be laid upon eT2{t) 
if it is shifted rightwards properly, i.e., the two curves 
are related as 



e[T2,Ti](i;iwi) ^ eT2{t2 + twf^) at t2 > t. 



/off') 
wl ' 



(3.2) 



with the same t^f ■* as the one determined before through 
i?[T2,Ti] (i; iwi)- Notice that the shifted e[T2,Ti](i; ^wi) 
approaches 67^2 {t) from above. 

In the positive T-s hift process, ejTTj.Ti] (^5 ^wi) behaves 
similarly and eq.(3^) also holds as shown in Fig. 4(b). 
In this case the original (shifted) e^j^^ j-^] (i; f^i) merges 
into (i) from above (below) . 

From inspection of Figs. 3 and 4 wc may divide the 
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Fig. 3. (a) Growth of ^2[T2,Ti](*; *wl) in a negative T-shift pro- 
cess with Ti = 0.8, T2 = 0.6 and twl = 1000 MCS (solid cir- 
cles) . The solid and broken lines are the isothermal aging curves 
Rxiit) and iJ^aCf), respectively. The open triangles represent 
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[T2,Ti] 



and t 



(cff) 



{t; twi ) plotted against t = t2 + t^^f ^ where t2 = t — t 



6230. (b) i?[T2,Ti] iwl) in a positive T-shift pro- 
cess with Ti = 0.6, T2 = 0.8 and iwi = 10000 MCS. The symbols 
and curves are the same as those in (a) but with t 
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aging process after the T-shift into the following three 
regimes. 

1) Quasi-reversal regime (ia < ti-cv — 1) 

Each spin responds individually to changes in its lo- 
cal Boltzmann weights. This gives rise to a relatively 
large change in e[T2.Ti]{t',twi) in the first MC step 
after the T-shift. The process is expected almost 
reversible, though we have not explicitly confirmed 
it. 

2) Transient regime (frcv ^ h ^ ^wi^) 

Judging from the behaviour of R[X2,Ti]it,twi), 
eq.(3.1), we may regard that the system has aged 



roughly for t^f'' at T = T2 already at around 
t2 ^ trcv Actually, ejT-^.TijC^; *wi) is significantly 
smaller (larger) than eT'2 (t) at i > iwi in the nega- 
tive (positive) T-shift process. 

Another interesting comparison is between 
e[T2.Ti]it',twi) as a function of t2 and eT2{t — twi"^) 
which is shown in the insets of Fig. 4. Now the for- 
mer is larger (smaller) than the latter in the negative 
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Fig. 4. Time evolution of the energy per density, e[T2,Ti] *wi )i 
in the same (a) negative and (b) positive T-shift protocols as 
in Fig. 3. The symbols and curves represent ejx-j.Ti] (*i *wl) 3,nd 
exit) in replacement of ijjy^ ,Ti] *wl) and RT{t), respectively. 
In the insets of both figures e[T2,Ti] *wi ) as a function of t2 
and e^jCt — ^wi^') shown by the data points and the curve, 
respectively. 



(positive) T-shift process, and it relaxes to the lat- 
ter slowly, by no means exponentially. Within our 
scenario introduced in §2.3, these can be explained 
qualitatively as follows. Just after the negative 
(positive) T-shift, the population of active droplets 
within the domain is larger (smaller) than in equi- 
librium at T2. Thus the energy density decreases 
(increases) toward the isothermal curves. 
3) Isothermal regime (^2 ^ ^wf^) 

In this regime time evolution of the system is noth- 
ing but isothermal aging at T = T2. The re- 
sult eg. ( p. 2D combined with eq.(3.1) means that the 
mean domain size and the energy density are related 
by the scaling relation eq.( |2.6D . The latter implies 
that relaxation of the energy density in this regime 
is due to the coarsening of domain walls. 
The existence of isothermal regime 3), to which the 



by the present simulation. It is certainly a consequence 
of the persistence of domains through the T-shift process 
described before, and can be regarded as one of the mem- 
ory effects. An interesting question is if the phenomena 
in transient regime 2) are related with the rejuvenation, 
or chaos effect observed experimentally. 

3.4 Spin auto- correlation function 

We define the spin auto-correlation function in the T- 
shift process as 



^[T2,Ti]{''"jtw2,'twl) — {Si{T + tyji + tw2)»S'i(^wl + ^w2))- 

(3.3) 

and measure relaxation of 

X[T2,Ti](T;i;iwi) = ^[1 - C[T2,Ti]{r;t2,t^i)] (3.4) 

with increasing t2 after the T-shift. At each measure- 
ment, a fixed value of the time separation r is chosen as 
a parameter. 

We naively expect that the FDT holds for time regime 
T < iw2 as in the case of isothermal aging (see eq.(2.11)). 
If it is the case the above X[T2,Ti] {t] t] twi) is identical to 
the ZFC susceptibility. In experiments, relaxation of the 
out-of-phase component of the ac susceptibility with a 
fixed frequency, say uj, is most frequently measured also 
in the T-shift protocol. It has been numerically ascer- 
tained in II that the relaxation of the ac susceptibility 
at a fixed frequency uj is very similar to that of the ZFC 
susceptibility with a fixed r = 27r/w = t^. Thus, in 
the present work, we regard that XlT2,Ti]i'''ui',t, twi) is es- 
sentially equivalent to the ac susceptibility measured in 
experiments. 

The function with t„2 ~ 0, denoted.by C[T2,Ti] {t; 0, iwi), 
was examined in our previous workES hereafter referred 
to as III. We found that at r ^ t'^f'^ it follows sim- 
ilar scaling forms to those for CT(T;iw) in isothermal 
aging described in §2.1 if the factor Lt{t) / Rritvj) in 
eq.( p.lOD is replaced by T[T2,Ti] ('r; iwi)/-RTi (iwi)- Here 
L[j-2 (r; iwi) is a mean size of the droplet excitations, 
or quasi-domains, which are in local equilibrium at the 
new temperature T2 within time scale r after the T-shift. 

In Fig. 5 we show typical results of X[T2,Ti]{Tuj]t',t^i) 
simulated for a fixed (= 32). Their behaviour is rather 
similar to that of ej^^.Ti] ^wi) shown in Fig. 4. Just 
after the T-shift, it rapidly decreases (increases) and then 
relaxes slowly to the isothermal curve XT2('^t»j;^) from 
below (above) for AT < (> ). F urthermore, it also 
satisfies a relation similar to eq.( |3.2| ): when the part of 
its curve at i2 > is plotted against i = t2+ i^f ■* it now 
relaxes to the isothermal curve from above (below). 



X[T2,Ti]{'ruj\t\t^i) 



XT2{Tu;t2 



at t2>t 



(off) 



wl 

(3.5) 

(~ 824) has been chosen so as to 
obtain a good collapse of the two curves at the time range 



Here, a value of i^f 



to > t 



(off) 



The resultant ■* is again consistent with 
combined with eqs.(1.2) and (1.3) obtained in 



eq.(p.l2| 
I. 



system crossovers from regime 2) at around ^2 ~ ^wf^ is feature in isothermal regime 3), ^2 > t 



(off) 
wl ' 



another important aspect of the aging dynamics observed 



means that there the mean domain size and the cor- 
relation function are related by the scaling relations 
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of eqs.(^), and ( |2.10D . Thus relaxation of 

X[T2,Ti] {tuj] t; twi) in this regime is due to increase of the 
effective stiffness due to the coarsening of domain walls. 

(a) 




Fig. 5. X[T2,Til {''"ij; *wi) in a negative T-shift with Ti = 
0.8, T2 = 0.6 (a) and in a positive T-shift with Ti = 0.6, T2 = 
0.8 (b). The symbols and curves, including those in the in- 
sets, represent X[t2,Ti] (''"ij ; *wl ) and XT{Ti^\'t) in replacement 
of e[T2,Ti] *wl) and ej'{t) in Fig. 4), respectively. The values 
of twi and t^f'^ are denoted in each figure. 



Behaviour of X[T2,Ti] ('^w; ^5 ^wi) in transient regime 2) 
is considered to reflect nature associated with quasi- 
domains. Let us here focus on the negative T-shift shown 
in Fig. 5(a). As pointed out just above, X[T2,Ti] (■^wi ^5 ^wi) 
is definitely smaller than XT2('''i^;0 a-t ^ ^ ^wi- This is 
different from the rejuvenation, or chaos effect observed 
in experiments if the latter means that the system after 
the T-shift looks almost as young as the system which 
has just started a new isothermal aging at T = T2. If 
this is the case, X[T2,Ti] [tuj', t', twi) as a function of t with 
t^f smallepJjhaTi t^i should coincide with {tlj ! t) 
at t2 > ti-cvE5> The behaviour of xji-^^j-j] (r^^; twi) may 
sug gest that, if the overlap length /at introduced by 
eq.( 2.13 ) in §2.3 exists in the present SG model, it should 
be larger than i?Ti(^wi) observed in the simulation. 

Let us then compare X[T2,Ti] {ti^', t; iwi) as a function of 



the latter and the excessive part decays very slowly. This 
feature in transient regime 2) may be explained within 
the scenario sketched in §2.3 as the following. After 
tw2 from the T-shift, there remains excessive population 
of droplets of length scales larger than the size of the 
quasi-domain Tjj-^ j^^] (^2, ^wi)- The excessive population 
is proportional to AT = T2 — Ti. The boundary of the 
latter droplets may act as frozen-in domain walls and 
reduces the effective stiffness of droplets smaller than 
L[T2.Ti]it2,tv/i) just as in the case of .the domain walls 
which separates different pure states.t3EJ) This yields 
the extra contribution, which is proportional to AT, to 
the relaxation. The scaling analysis on the excessive part 
in terms of T[T2,Ti] {h, twi) and i?[T2,Ti] {t'-, twi) is now un- 
der investigation. 

3.5 T -cycling process 

A T-cycling process is a combination of two T-shift 
processes with the same |AT| but with opposite signs. 
At a time ^2 = ^w2 after the T-shift from Ti to T2 the 
temperature is turned back to Ti, and the subsequent 
relaxation is observed. When t^2 ^ ^wf^ the T-change 
back to Ti is effectively the same as a simple T-shift 
process from an isothermally aged state at T2 to Ti, 
as shown in Fig. 6(a). In this negative T-cycling pro- 
cess (T2 < Ti), the approximated ac susceptibility at 
t = t3 + t^2 + tvji, now denoted simply by XT{t){'''uj',t), 
merges into the isothermal curve XTi{Tco]t) from below 



(efl) 



Here 



t2 with XT2 {tlu', ^"^wf )• The two curves are shown in the 
inset of Fig. 5(a). Interestingly, the former is larger than 



when the former is plotted against t = t:} + t 
tw*'' is the effective waiting time given by £^2^^^'^^'^'^^ 
(since iw2 3> i^f"* in this case). In the case that tw2 
is small enough, we obtain a similar behaviour but with 
tw^^ — iwi (not shown). This means that aging at T2 
does not contribute at all to aging at Ti . This is nothing 
but the mernp|K Jiffect and has been observed by many 
experiments.o' LJ'Es' 

In a positive T-cycling process with T2 > Ti shown 
in Fig. 6(b), our simulated data still exhibit the mem- 
ory effect in the sense that the curve at ^3 > lies on 
the isothermal curve of Ti if it is shifted by a proper 
amount, i.e., t^w^^ which depends not only iw2 but also 
on tvii- This type of memory retained over a positive 
T-cycling has been ipdeed observed in the ferromagnetic 
fine particle system. Ej-* 

3.6 Continuous T-change with an intermittent stop 

We have also performed simulations on a continuous 
T-change with an iiitermittent stop. Following the ex- 
perimental protocolj23) we cool a system by a constant 
rate from above Tc to below with an intermittent stop 
at a certain temperature Ti below Tc, restart the cool- 
ing, and then heat it back by the same rate as in the 
cooling. A typical result is shown in Fig. 7. In this 
analysis XT{t){'^ui]t) has been directly measured as the 
response to an ac field h{t) = Hq sin(wt) with Hq = 0.3 
and 27r/w = 160 (MCS). 

At the intermittent stop at Tj, x'T(t)i'^'^^) relaxes 
downward and deviates from the reference curve of cool- 
ing. After restarting the cooling, XT(t)('^j*) i^ seen to 
merge into the reference curve. It appears that the long 



(a) 



Tatsuo KOMORI, Hajime YOSHINO and Hajime Takayama 

[xlO-'] 





0.3 




0.28 




0.26 












0.24 










0.22 








n 9 
u.z 




0.18 




0.16 




0.3 




0.28 




0.26 












0.24 










0.22 








0.2 




0.18 




0.16 


Fig. 


6. X 


0.8 


, T2 


0.6 


, T2 





?wi=512, rw2=16384 
;r«=2411 




\^ ri=o.8 


r,=o.6\ 




Tb=27i/g)=32 





100 

(b) 



1000 ~ 



10000 



100000 



?wi=512, rw2=512 
?r'«=5792 




ri=o.6 \ 



Tb=27i/g)=32 



100 1000 ~ 10000 



100000 



(off) 



in Fig. 5 but with t = ts + t 



time aging at Tj does not affect subsequent cooling at 
temperatures lower than Tj by a certain amount. Here 
we call this phenomenon a rejuvenation-like effect, since 
X^^^^(u;;t) merges with the reference curve but without 
a strong vertical increaseJji contrast to the rejuvenation 
seen in the experiment .tS' By a detailed inspection of 
Fig. 7, one may notice that the curve XT{t)('^'T^) looks 
as if it increases just after the cooling is restarted. We 
consider that this behaviour does not indicate the reju- 
venation but is due to the following reason: ^^^^^(a;;^) 
during the stop initially decreases appreciably but finally 
its relaxation becomes invisible because of its statistical 
fluctuation. The curve looks as if it increases initially 
at the restart of cooling simply because the relaxation 
starts in the middle of the fluctuation. 

When the system is heated back from the lowest tem- 
perature of cooling process, x'T{t)i^'^) exhibits a dip at 
around T — Ti, though it is rather shallow. This means 
that the memory imprinted by the long time aging at 
Ti in the cooling process, in fact, has not been erased 
until the system comes back to Ti from the lower tem- 
peratures. To conclude the present SG model exhibits 




Fig. 7. XT(t) {''''^ ; in a continuous T-change with an intermit- 
tent stop. We cool and re-heat the system continuously with 
a constant rate of J/3200 per MCS from T = 3.0 to T = 0.1. 
The open circles represent the reference curves in cooling (upper 
curve) and re-heating (lower curve) protocols. The solid circles 
represent the same protocol except for an intermittent stop in 
cooling at Ti = 0.8 for 32000 MCS. 



the rejuvenation-like and memory effects within a time- 
window of the present simulation. 

§4. Discussions 

We have studied aging phenomena in various T- 
change protocols on the 3D Gaussian Ising SG model 
by Monte Carlo simulations. In particular, we have ar- 
gued that, in a T-shift protocol, slow aging dynamics 
after the T-change from Ti to T2 at ^2 = consists 



transient {t2 < t)^i ) and isothermal 



(^2 ^ ^wf ) regimes, where t^f'' is the effective wait 



of two regimes: 

> 

ing time specified by eq.( ^.12 ). In the transient regime, 
re-distribution of thermal weights, from those at Ti to 
T2, of droplet excitations with the characteristic length 
scale L[T2,Ti]{t2, ^wi) is taking place within each domain 
of the mean size i?[T2,Ti] ^wi)- We called the length 
scale i[T2,Ti] (^2, ^wi) as the size of quasi-domains. 

Suppose, in a negative T-shift process, we further de- 
crease temperature to T3 (< T2) at t2 = tw2 *wf^) 
and wait a certain period. Then there appear now quasi- 
domains locally equilibrated to T3 within each quasi- 
domain of T2. Because of the wide separation of time 
scales at different length scales emphasized in §2.1, a 
nest (or hierarchy) of quasi-domains can be realized dy- 
namically when the system is cooled step-wise, or con- 
tinuously. At low temperatures the ac field can only 
excite droplets of the size Lt{1/oj) corresponding to the 
frequency lu. If the size is small enough, a new contri- 
bution to the ac response will appear every time a new 
frozen-in domain wall appears associated with the quasi- 
domain due to a step of cooling. This picture, quasi- 
domains within (quasi-) domains, is our interpretation 
of the rejuvenation-like effect seen in Fig. 7. Thi? picture 
is somewhat similar to the one proposed before.E3'E3-' By 
our interpretation the observation of the rejuvenation- 
like effect crucially depends on time and length scales by 
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which one looks at the system. It is the easier to be ob- 
served, the shorter and smaller are the time and length 
scales. 

By means of our scenario the memory effect is rather 
easy to be understood. For example, the continuous T- 
change protocol shown in Fig. 7 is interpreted as follows. 
When the system is re-heated from the lowest tempera- 
ture, the hierarchical structure of quasi-domains is erased 
from the lower levels. When the temperature comes back 
to Ti, the large domains imprinted at the intermittent 
stop in the cooling process reappear and govern the re- 
sponse of the system, which causes the memory effect 
in this protocol. According to our scenario, a memory 
imprinted by aging at a certain temperature Ti is erased 
when the system is aged at different temperature, say 
T(j, for such a long period that the mean size of quasi- 
domains of Tfi catches up the domains imprinted at Ti, 
as what happens in regime 3) of the T-shift protocol. 
Even in that case, the memory of the previous thermal 
history is kept in the effective waiting time, i.e., the sys- 
tem looks more aged than the period that it is actually 
aged at T = Td (see eq.(U)). 

Let us now compare our simulated results with exper- 
imental observations. As pointed out already in §1, in 
the ferromagnetic fine particles system (FFPS), whose 
aging dynamics has been recently studied by Mamiya et 
al.,cB the microscopic time, Tj^ic, required for a mag- 
netic momeat of each fine particle to flip is of the order 
of 10~'^sec.E3' Therefore the time window of their obser- 
vation, 1 ~ 10"^ ksec, in unit of Tmic is not much different 
from that of our MC simulation. It is also noted that 
in the FFPS each magnetic moment behaves as an Ising 
spin due to strong magnetic anisotropy. 

Indeed in the T-shift protocol, Mamiya et al. have 
observed relaxation of X[T2 Ti] ("^j ^' ^wi) which is quali- 
tatively quite similar to our simulated X[T2,T{\{Tu]t',iwi) 
shown in Fig. 5. They have determined amount of the 



peak position of its relaxation rate with respect to Int 
may be assigned as i^f ■* (crossover from regime 2) to 
regime 3)). The experimentally observed i^f'* becomes 
larger (smaller) than t^i in the positive (negative) T- 
shift, which is qualitatively consistent with eq.( ^.12 ) with 
monotonic increase of l/z{T) on T. Quantitatively, how- 
ever, the shift observed experimentally is much larger 
than what is expected from eq.(1.3). 

Djurberg et alBS> furthermore claimed that they ob- 
served the chaos effect in T-shift processes with |AT| 
larger than a certain value. The rejuvenation, or chaos 
effect in the T-shift and T-cycling protocols obsjetv 
in XT(t)('^'^) been reported for CdCri.7lno.3S4ll3'l 
and Cu:Mn 2at%ii' spin glasses, as already—mentioned 
in §1. The observed XT{t)(^'^^) claimedEZP to indi- 
cate that the system is literally rejuvenated in the sense 
described in §3.4. Neither such a chaos effect nor lit- 
eral rejuvenation have beep_seen in our present simula- 
tions (see also Kisker et alS3 ) . At the moment it is not 
clear yet whether this discrepancy originates simply from 
the difference in time scales of the observations, or from 
something else. 

To conclude, we have simulated various T-change pro- 
tocols of aging in the 3D Ising SG model, and have 
proposed a scenario in terms of domains and quasi- 
domains: the domains ever grow in any aging process 
in the SG phase and the quasi-domains appear within 
the domains or within themselves. The results observed 
in our simulations are reasonably well interpreted by the 
scenario. Qualitatively, it is also consistent with many 
features observed experimentally, particularly, those re- 
lated to the memory effect. For certain phenomena such 
as the growth law of Rrit), the consistency is even semi- 
quantitative. Concerned with the rejuvenation, or chaos 
effect, however, the scenario may not be applicable, and 
further studies are required. 



shift of x(t2,Ti](^;^;*wi), 



or t^f\ by the same method Acknov^fledgments 



that we have used for analysis of X[T2,Ti]{Tuj',t;tv,i) in 
Fig. 5. The effective waiting time they obtained is 



written as t 



trast, t^^^^ extracted from our simulation is given 



(off) 
wl 



3 X (iwi)'' with b ~ 1.0. 



By con- 

by 

eg. ( 2. 12 ), which, with making use of eq.(1.2), is reduced 
to t^(x (twi)'' with b = z{T2)/z{Ti) (~ 1.11 (1.25) for 
Ti = 0.8 and T2 = 0.7 (0.6)). Also aging dynamics of 
the FFPS in the T-cycling protocol agrees qualitatively 
with the results of our simulation. In particular, the ef- 
fective waiting time iw^'' which depends on t^i (memory 
effect) is observed even in the positive T-cycling process, 
though Mamiya et al. have claimed that the memory is 
partially lost because some particles with a larger size 
than the average are affected by the chaos effect associ- 
ated with the temperature change. 



For ordinary spin glasses with Tmic ~ 10~ "^sec, there 
have been only a few experiments, in which t^^^^ in thCf^j- 
shift process is systematically studied. Djurberg et al& 
measured the ZFC magnetization of Cu:Mn 13.5at% spin 
glass with a fixed t^i and varying AT. Via the FDT 
similar to eq.(2.11) their analysis is directly compara- 
ble with our simulation on Cittj.Ti] (''"; 0; ^wi) h^ HI. The 
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By means of Monte Carlo simulations on the three-dimensional Ising spin-glass model, we have 
studied aging phenomena with various temperature(T)-change protocols. Particularly, a T-shift 
protocol, in which a system is first quenched to and aged for a period twi at a temperature Ti, 
and subsequently aged at a new temperature T2 is closely investigated. Most importantly, the 
mean size of domains, which is extracted from the replica-overlap function, was found to grow 
monotonically without any appreciable decrease by the T-change. We also found the relaxation 
of energy density and spin auto-correlation function during the T-shift process can be explained 
within the following picture: the dynamics finally crossovers to an isothermal aging of T2 at 
around t2 — t^f after the T-change, where t^f is an effective waiting time related with twi, 
but in the transient regime t2 ^ ^wf i adjustment of the population of thermally active droplets 
from that at Ti to T2 takes place very slowly. Implications of the results on aging phenomena in 
other T-change protocols such as T-cycIing and continuous T-change with an intermittent stop 
observed by simulations as well as experiments are discussed. 
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§1. Introduction 

Aging phenomena in spin glasses have been extensively 
studied in recent years. Experimentally, many elabo- 
rated protocols have been..adopted to reveal various as- 
pects of the phenomena.El'13'El' Theoretically, on the other 
hand, many ideas and models on spin glasses and re- 
lated systems have been proposed and examined to get 
a proper_|iiiderstanding of aging dynamics in such glassy 
systems .□'B-' On the side of numerical study, a number of 
simulations, which could account for quali1(atii2.e features 
of the aging effects, have been performed .oQ'Ij'' In spite 
of the enormous efforts, however, a unified picture of ag- 
ing dynamics in spin glasses has not been established 
yet. The difficulty is directly connected to nature of the 
low-temperature spin-glass (SG) phase which has been a 
controversial issue since more than a decade ago. 

Among simulational results on the short-ranged Ising 
SG models, of particular importance is the growth of a 
mean domain size which is observed through the replica- 
overlap function, G{r,i), in isothermal aging after a sys- 
tem is quenched instantaneously from temperatiii:e.abj:jaa3 
its SG transition temperature to T below Tc.ollj'lllltH' 
The function G{r,i) is defined as 

i = l ™ 

where a and [3 are indices of two replicas which have dif- 
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ferent random spin configurations at t = (an instant of 
the quench), and are independently updated. Its correla- 
tion length Rxit) is found to grow almost ma_nQwer-law 
in the case of 3 dimensional (3D) models,t3Eil'E3'' 

Rrit) - Lo{t/Toy/'^^\ (1.2) 

where Lq and tq are certain characteristic units of length 
and time, respectively, and the exponent l/z{T) depends 
on T. In our previous work,llj'' hereafter referred to as I, 
we obtained 

l/z(T) ~ 0.17T, (1.3) 

at T < 0.7J for the Gaussian model with zero mean 
and variance of interactions J. More precisely, it can 
be also fitted to a logarithmic functional form written 
as i?T(i)_£i_/?o + b(\niY^^ with Ro,b and tp being con- 
stants.BEJ'll^) This ambiguity is attributed to the insuf- 
ficiency of the time window of simulations by presently 
available computers. Although there remains such an 
ambiguity in the explicit growth form of Rxit), we can 
naturally regard Rxit) as a characteristic length scale of 
the SG ordering in aging dynamics, i.e., a mean distance 
of domain walls separating different pure states. 

It is natural to ask to what extent the mean domain 
size Rxit) is related to aging behaviours observed in-yaj:- 
ious quantities simulated. The droplet theory,ll3'Mr^lEZP 
particularly the one proposed by Fisher and HusaiZI* con- 
tains scaling ansatz which relate various quantities to 
the characteristic length scales such as Rxit) mentioned 
above. We have been numerically studying aging phe- 
nomena in thp_3Jl Ising SG model to clarify the validity 
of the ansatz.BBy) 

Experimentally, various interesting phenomena have 
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Fig. 1. Schedules of T-change in some experimental protocols of 
aging: 1) isothermal, 2) T-shift, 3) T-cycle, and 4) continuous 
T-change with an intermittent stop. 



been also observed in aging processes other than the 
isothermal one. They are T-shift, T-cycling, and con- 
tinuous T-change with an intermittent stop(s) protocols 
whose schedules of T-change are schematically shown in 
Fig. 1. The ac susceptibility has been continuously mea- 
sured through these T-changes. In the case of a negative 
T-shift process with AT = T2 — Ti < in Fig. 1 and 
with a magnitude |AT| larger than a certain value, the 
ac susceptibility is observed to increase discaotiDUDusly 
just after t = t^i and then to relax again. tSEJ'EiJ) Its 
behaviour looks as if the system restarts aging by the 
T-shift at t = tvfi forgetting its thermal history before 
twi- It is called the rejuvenation, or chaos effect. Fur- 
thermore when the temperature is turned back to Ti at 
t = tmi+^w2 (T-cycling protocol), the value of the ac sus- 
ceptibility returns to that on the isothermal aging curve 
of Ti at almost t = twi, and then relaxes along the curve. 
This behaviour indicates that the system in fact remem- 
bers its history before twi, and is called the memory 
effect. Similar coexistence of the paradoxical phenom- 
ena, rejuvenation and memory effect, has been recently 
observed also in— a-continuous T-change with an inter- 
mittent stop(s).E3E3' These ejcperimental findings have 
stimulated theoretical interest .eP 

The main purpose of the present paper is to investi- 
gate aging phenomena in various T-change protocols by 
simulations on the 3D Gaussian Ising SG modelEil* We 
concentrated mostly on the T-shift process, which is the 
most fundamental among the various T-change proto- 
cols, since the latter are considered as certain combina- 
tions of the former. We investigated the time evolution 
of the mean domain size and the relaxation of the energy 
density and the spin auto-correlation function during the 
T-shift process. 

We have found, within the time window of our simu- 
lations, the following characteristic features of the aging 
dynamics, i) The mean domain size R[T2,Ti]{i] ^wi) in the 
T-shift process (see Fig. 1), with either a negative or pos- 
itive AT, does continue to grow. Here and hereafter t is 
a time measured from the first quench from above Tc to 
Ti, t'z = t — t^i a period that the system ages at T = T2 
after the T-shift at t = t^i . ii) The time evolution of 
the energy density and the spin auto-correlation func- 



tion is related with the mean domain size -R[T2,Ti](^; ^wi) 
at sufficiently long time scales after the T-shift by the 
scaling relations which are the same as isothermal aging, 
iii) However, there is a transient time regime after the 
T-shift where there is extra contribution to the relax- 
ations. Interestingly, the extra contribution also shows 
very slow, non-exponential relaxation. We interpret the 
latter as not due- .to the chaos effect predicted by the 
droplet theorytStZP but as due to adjustment of the pop- 
ulation of droplet excitations within each domain. 

The memory effect observed experimentally is easily 
interpreted from the above-mentioned scenario, namely, 
it is attributed to the persistence of domains in a T-shift 
process. It is not yet conclusive, however, whether na- 
ture of the transient regime in the scenario is common to 
phenomena observed experimentally just after a T-shift. 
In this respect, the recent work on aging in the ferro- 
magn£lic fine particles system (FFPS) done by Mamiya 
et alx3 is of quite interest, since the time window of their 
observation in rUJiit of the microscopic time (~ 10~^sec 
for this systemE^) is rather closer to that of our simula- 
tions. Indeed, they have observed T-shift aging phenom- 
ena which can be well interpreted by our scenario. It 
seems, however, that the scenario is hard to explain the 
rejuvenation, or chaos effect. observed experimentally in 
ordinary spin glasses,ll3'Eilt3E3' whose microscopic time 
is of the order of 10~^^sec. 

The present paper is organized as follows. In the next 
section we briefiy review the scaling properties of isother- 
mal aging derived from the droplet theory, add further 
comments on them, and present the scenario on the T- 
shift process. In §3 we present the results of our simula- 
tions on various aspects of the T-shift, as well as the T- 
cycling and continuous T-change processes. In the final 
section implications of the results with the experimental 
observations are discussed. 

§2. Phenomenological Picture 

2.1 Scaling properties of isothermal aging 

Here we briefiy review the scaling propertie.s nf isother- 
mal aging derived by the droplet theoryllJ'li3 due to 
Fisher and Huse (FH). In an isothermal aging process 
a system is quenched from above Tc to temperature T 
below Tc. At waiting time t^ after the quench, there 
are domain walls separating different pure states of the 
SG phase at a typical distance Rxitvi) from each other. 
The distance Rxitvi) can be measured in simulations for 
instance through the replica overlap function defined in 
eq.(O). 

Except for places where the domain walls run, the bulk 
of the system is essentially in equilibrium. An impor- 
tant subtle feature inside the domains of a spin glass is 
that there can be large scale thermal fiuctuations due to 
droplet excitations. In ideal equilibrium, a droplet exci- 
tation can be defined as a global fiip from a ground state 
of a droplet (cluster) oLapins within a distance L/2 from 
a certain given site i.t^P The latter can be considered 
as a simple two-state system with a free-energy excita- 
tion gap Ti(i) and a barrier energy 5l(«) of a thermal 
activation process. 
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The typical value of the gap scales as, 



ptyp 

^ L 



r{L/Lo) 



(2.1) 



where T is the stiffness constant of domain walls on the 
boundary of the droplet. The gap is however broadly 
distributed and its distribution function Pl{F) is con- 
sidered to follow the scaling form, 



Pl{F) 



1 



F 



(2.2) 



A very important property of the scaling function p{x) 
is that it has finite intensity at x = 0, p{0) > 0, which 
allows large scale droplet excitations even at very low 
temperatures. The probability to have a thermally active 
droplet at a very low temperature T is given by 



p{L;T) 



rP(0). 



(2.3) 



The latter feature yields equilibrium properties which 
make spip_|glasses distinctly different from simple ferro- 
magnets.tS 

The relaxation time of a droplet of size L centered at 
site i is given by the Arrhenius law 

'BL{i) 



To exp 



(2.4) 



where B^li) is the energy barrier of the droplet. The 
typical value 5^^'^ of the energy barrier B^li) scales as. 



Ejtyp 



A{L/Lo 



(2.5) 



where A is a characteristic unit for energy barriers as- 
sociated with thermal activation of droplet excitations, 
and the exponent ip satisfies 6 < ip < d — I. 

In the droplet theory, it is assumed that the growth 
of domains is also governed by droplet excitations. An 
obvi ous but very important consequence of eqs.(2.4) and 
(2.5) is that dynamical processes, including both the do- 
main growth and droplet excitations, ajt different length 
scales are extremely widely separated.EP Hence within a 
time scale of small scale processes, large scale processes 
look as if they are almost frozen. 

In isothermal aging, for example, the relaxation of the 
excessive energy per spin Sexit) with respect to the equi- 
librium value is expected to follow the scaling form as. 



Serit) ~ T{RT{t)/Loy/{RT{t)/Lo) 



(2.6) 



where T is a characteristic energy scale, and d is the 
dimension of space which is 3 here. This scaling form 
was confirmed, in our previous work I, with the energy 
exponent 6 = 0.20 ± 0.03 which agrees with the result of 
the defect energy analysis at T = OE£P 

The droplet theory also predicts important scaling 
properties of the spin auto-correlation function. 



Ct(t;K) = {Si(T + K)Si(t^)), 



(2.7) 



in isothermal aging. Here the over-line denotes the av- 
erages over sites and different realizations of interactions 
(samples) , and the bracket denotes the average over ther- 
mal noises. In the so-called quasi-equilibrium regime de- 
fined by T <^ tvf, the typical size Lt{t) of droplet exci- 



tations which take place in the time scale of r is much 
smaller than the typical separation i?T(tw) of domain 
walls present after the waiting time For this situ- 
ation the droplet theory introduces a phenomenological 
concept called effective stiffness which characterizes a re- 
duction of the excitation gap of small scale droplets due 
to the presence of the domain walls. The latter yields 
tw-dependence of CT(T;tw) which is written in terms of 
-Rt(^w) and Lt{t) as the following: 



with 



1 - Ct (r ; tw) = ACt (r ; t^) + ar (r) , 



arir) = 1 - C^^^\r) 



(2.8) 
(2.9) 



and 



T 



Lo \ f Lt{t) 



d-0 



'T \Lt{t)) \RT{t^), 

(2.10) 

Then the zero field cooled susceptibility XzfC;t(''"; ^w) 
which is often observed by experiments should also be 
understood within the same ansatz since the fiuctuation- 
dissipation theorem (FDT) 

XzFC;T(r;tw) ~ ^[1 - CT{T;t^)] = XT(r;tw). (2.11) 

holds in the quasi-equilibrium regime. The FDT 
has bean— tested and confirmed by the careful experi- 
menti£3'LlP even in the presence of weak violation of the 
time translational invariance. 

2.2 Further comments oil isothermal aging 

In our previous workJ13 hereafter referred to as II, 
we confirmed the scaling ansatz on CTiT'jt^) described 
in the previous subsection. Here, we first make the com- 
ment that XzFC:T (t; ^w) experimentally observed and an- 



alyzed through eq.( 2.11 ) also follow the scaling ansatz. 

In Fig. 2, for example, we plot XzfC;T (''"; ^w) of a 
Ag:Mn 2.6at% spin glass against t^/r. The latter is 
evaluated from niTRM;T (''"; ^w) measured by the Saclay 

grouplil) as XZFC;T(T;tw) = [niFC;T - niTRM;T(T;tw)]//j 

with mYc,T being the field-cooled magnetization under a 
field of strength h. The data for different r's lie on a uni- 
versal curve when they are vertically shifted by suitable 
amounts depending on r as shown in the inset of the fig- 
ure. The feature is quite similar to 1 — CT(T;tw) analyzed 
in II. [Note that the factor (Lq/ Lt(t)Y is practically 
constant in eq.( 2.10| ) since 9 is much smaller than d — 0, 
as was the case of 1 — C t{t; tw) in H.] This indicates that 
XzfC;t(''"; ^w), via eq.(2.11), follows the scaling forms 



eqs.( |2.8| ), ( |2.9| ) and ( 2.10| ) as well. The values of the expo- 
nents {d — 6)/z{T) and 6/z{T) turn out to agree roughly 
with the estimates using 6 and_l/z(T) ~ 0.18(T/Tc), i.e., 
eq.( p_.3| ) with Tc = 0.95±0.04j23' extracted from our anal- 
ysis in II. However, such semi-quantitative agreement is 
lost for the data at temperatures closer to Tc. 

An intriguing point is that the above result sug- 
gests that Rxit) in the time scale of experiments is 
also compatible with the estimate from our simulations. 
In_Lhis context it is worth pointing out that Joh et 
alS3^ have recently extracted Rxit) of Cu:Mn 6at% and 
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Fig. 2. xzFC;T(^;tw) of a AgMn spin glassEf at T/Tc = 0.87 
for fixed r's plotted against tv^/r; r = f , 2, 4, 8, 256sec from 
bottom to top and tw = 300, 1000, 3000, 10000, SOOOOsec. In the 
inse t the scaling behaviour of AxzfC;t('^5^w) corresponding to 
eq.( ^.lo| ) is demonstrated. 



CdCii 71110.384 spin glasses from the thermoremanent 
magnetization (TRM). The growth law they have ob- 
tained turns out to agree with the simulational result of 
eqs.(1.2) and (1.3) even quantitatively. The estimated 
Rxit) with t of the order of laboratory time is at most 
only several tens of lattice distances. However, the ex- 
perimental studies on the isothermal aging dynamics so 
far reported have not converged to a unified picture yet 
as we discussed in II. 

Another comment is related to the characteristic 



length and time scales, Lq and tq in eq.(L2), which we 
did pjat fix explicitly in II. According to the droplet the- 
ory|13 they are some microscopic scales at lower temper- 
atures, but near below Tc they are expected to crossover 
to the SG coherence length ^_ and the corresponding 
criticairelaxation time, respectively. Indeed, Hukushima 
et a/.ES' have found recently such a crossover behaviour 
in aging dynamics in the 4D ± J Ising SG model. A sim- 
ilar analysis on the present 3D model certainly has to be 
done, which may clarify the disagreement between the 
experimental and simulational results close to Tc jaaen- 
tioned just above. Furthermore, as suggested in,Ej) it 
may also provide a useful key to resolve the apparent 
ambiguity of the growth law mentioned previously. 

2.3 T -shift process mthm overlap-length 

Now let us consider the T-shift process. After a system 
is first quenched from above Tc to a temperature Ti be- 
low Tc, domains of different pure states grow up. Their 
mean size reaches to -Rti(^wi) after a given waiting time 
t = tvfi. After the temperature is changed at t = t^i to 
a new temperature T2, we expect the domains continue 
to grow but with a rate specific to the new temperature 
T2 (see eq.(3.1) below). It can be said that, at the in- 



stant of the T-shift, the system is aged as if it is aged 
for the effective waiting time t[^f^ by isothermal aging 



is determined by 



(2.12) 



with Rxit) being the growth law of the mean domain 
size in isothermal aging at T. 

Another important process which must take place after 
the T-shift is adjustment of the population of thermally 
active droplets within the domain. Since the equilibrium 
population is proportional to T as given in eq.( ^.3| ), there 
must be a change of the population of an amount pro- 
portional to AT = T2 — Ti at each length scale L. It 
is emphasized that this process itself is intrinsically very 
slow since adjustment of the droplets takes place only 
by activation processes described by eqs.(2^) and (3.5). 
Thus we naturally expect that there is a new charac- 
teristic length i[T2,Ti](^2,^wi) which slowly grows after 
the T-shift such that the population of droplets smaller 
than i[T2,Ti] (^2, ^wi) is equilibrated to the new tempera- 
ture T2. However, the population of droplets larger than 
-^[T2,Ti](^2,^wi) is still equilibrated with respect to the 
old temperature Ti . We call the characteristic length 
-^[T2,Ti](^2,^wi) as a size of quast-domams. 

The growth law of i[T2,Ti](^2, ^wi) is expected to be 
the same as the domain. When the system ages at T2 for 
a period about '' defined by eq.(2T2), i[T2,Ti](^2, ^wi) 
catches up -R[T2,Ti](^; ^wi) and the quasi-domains merge 
into domains in isothermal aging at T2. Thus at large 
time scales ^2 ^ ^wf'' i aging becomes essentially the same 
as isothermal aging at T2. The time range ^2 ^ ^wf"'' °^ 
the other hand, is regarded as a transient regime where 
adjustment of the population of droplets are taking place. 

In the above argument we implicitly assumed that 
equilibrium states at different temperatures are the same 
except for the change of the population of droplets 
eq.(2.3). It should be remarked that the assumption 
makes sense only for lenalii_spales smaller than the so- 
called overlap length /at-^SH' It is predicted that equi- 
librium states at different nearby temperatures Ti and 
T2 differ from each other at length scales larger than 
this length ^A'p-i This aspect is called chaotic nature of 
the SG phascL^ The overlap length /at is expected to 
diverge when AT = T2 — Ti — as 

2l(d,-20) 



Li 



VTi/2|AT| 



(2.13) 



the fractal 



where, 6 is the exponent in eq.(2.1) and 
dimension of surface of the droplets. 

§3. Results of Simulations 

3.1 Model and method 

We have carried out MC simulations on various aging 
processes in the same 3D Ising SG model as in our previ- 
ous works, i.e., the one with Gaussian nearest-neighbor 
interactions with mean_jiero and variance J = 1, and 
with Tc = 0.95 ± 0.04.y) The heat-bath MC method 
we use here is also the same as in I and II. The data we 
will discuss below are obtained mostly at T = 0.6 ~ 0.8 
in systems with linear size = 24 averaged over a few 
hundreds samples with one MC run for each sample. In 
I, it was confirmed that finite size effects do not appear 
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for these parameters within our time window ( ^ 2 x 10® 
MCS). 

3.2 SG coherence length — domain growth 

Let us begin with time evolution of -R[T2,Ti](^;^wi), a 
mean domain size in a negative T-shift protocol with 
AT = T2 — Ti < (see Fig. 1). The simulated results 
with Ti = 0.8, T2 = 0.6 and t^i = 1000 MCS are pre- 
sented by solid circles in Fig. 3(a). They are extracted 



from the replica overlap function eq.(l.l) which is ob- 
tained by averaging over several thousands samples with 
L, = 32. 

Clearly, -R[T2,Ti](^; ^wi) does grow continuously with- 
out any noticeable decrease by the T-shift. This im- 
plies that, just after the T-shift, the system looks as 
if it has aged isothermally at T = T2 for a period of 
the effective waiting time t^f'' defined in eq.(2.12). In- 
deed, when R[T2,Ti]{i]'tvii) at t > twi is plotted against 
i = t^f'' +t2 (open triangles in the figure), it almost lies 
on the isothermal aging curve Rx^it), i.e., the two curves 
are related as 



-R[T2,Ti](^;^v 

Here, a value of t 



(eff) 



at t2 ^ 0. (3.1) 
6230) has been chosen so as to 



obtain a good collapse of the two curves at time range 
i'j <^ ^w?- T^^^ chosen value turns out to satisfy eq.(2.12) 



combined with eqs.(1.2) and (1.3). 

As shown in Fig. 3(b), -R[To.Ti](^; ^wi) in a positive 
T-shift protocol is also described by eq.(3T), but with 
(- 1400) < twi (= lO'*) in this case. The contin- 
uous growth of domains in the T-shift process, either a 
negative or positive one, is one of the most important 
results observed in the present simulation. 

3.3 Evolution of energy — Transient and isothermal 
regimes 

In Fig. 4(a) we present the time evolution of energy per 
spin, e[T2_Tj](t; twi), in the same negative T-shift process 
as shown in Fig. 3(a). The original data drawn by the 
solid circles quickly deviate from the isothermal aging 
curve, eTj(t), just after the T-shift at t = twi, and then 
slowly crossover to ex^it) from below at longer times. 
As drawn by the open triangles in the figure, the part of 
e[T2,Ti]{i]'tvii) after the crossover can be laid upon ex^it) 
if it is shifted rightwards properly, i.e., the two curves 
are related as 



e[T2,Ti]{i;iwi 



at ^2 ^4f\ 



(3.2) 



with the same t[^f^ as the one determined before through 
-R[T2,Ti](^; ^wi)- Notice that the shifted e[T2,Ti]i't','twi) 
approaches eT2{t) from above. 

In the positive T-shift process, ej^^ ^j] (t; twi) behaves 
similarly and eq.(3.2) also holds as shown in Fig. 4(b). 
In this case the original (shifted) e[T2_Tj](t; twi) merges 
into eT2{t) from above (below). 

From inspection of Figs. 3 and 4 we may divide the 
aging process after the T-shift into the following three 
regimes. 

1) Quasi-reversal regime (^2 < ^rev — 1) 

Each spin responds individually to changes in its lo- 
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Fig. 3. (a) Growth of ^*^[T2,Ti] (^j ^wl ) in a negative T-shift pro- 
cess with Ti = 0.8, T2 = 0.6 and twi = 1000 MCS (solid cir- 
cles). The solid and broken lines are the isothermal aging curves 
RT-^{t) and Rj'^^t), respectively. The open triangles represent 
/?[T2_Tj](t;twl) plotted against i = t2 + where t2 = t — ty,i 



id t 



(eff) 



— 6230. (b) R[T2,Ti]{'^'i'^wi) in a positive T-shift pro- 
cess with Ti = 0.6, T2 =0.8andtwi = 10000 MCS. The symbols 



and curves are the same as those in (a) but with t 
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1400. 



cal Boltzmann weights. This gives rise to a relatively 
large change in ej^^ ^j] (t; twi) in the first MC step 
after the T-shift. The process is expected almost 
reversible, though we have not explicitly confirmed 
it. 

2) Transient regime (tj-ev <, t2 <, t[^f^) 

Judging from the behaviour of -R[T2,Ti](^; ^wi), 
eq.(3.1), we may regard that the system has aged 
at T = T2 already at around 



roughly for t^^f* 



t2 ~ tj-ey Actually, e[X2,Ti](i]'twi) is significantly 
smaller (larger) than ey^ (t) at t ^ twi in the nega- 
tive (positive) T-shift process. 

Another interesting comparison is between 
e[T2,Ti]{i','tvii) as a function of ^2 and eT2{t - ^^f) 
which is shown in the insets of Fig. 4. Now the 
former is larger (smaller) than the latter in the neg- 
ative (positive) T-shift process, and it relaxes to the 
latter slowly, at least not exponentially. Within our 
scenario introduced in §2.3, these can be explained 
qualitatively as follows. Just after the negative 
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Fig. 4. Time evolution of the energy per density, e[7^2,Ti] (^5 ^wl )? 
in the same (a) negative and (b) positive T-shift protocols as 
in Fig. 3. The symbols and curves represent e[7^2,Ti] (^j ^wl ) and 
e7^(t) in replacement of ^[T2,Ti] (^j ^wl) and Rj^^t), respectively. 
In the insets of both figures ^[t^^ ,Ti] (^j ^wl) as a function of t2 
and 67^2 ~ ^wi^^ ) shown by the data points and the curve, 
respectively. 



(positive) T-shift, the population of active droplets 
within the domain is larger (smaller) than in equi- 
librium at T2. Thus the energy density decreases 
(increases) toward the isothermal curves. 
3) Isothermal regime (^2 ^ ^w?) 

In this regime time evolution of the system is noth- 
ing but isothermal aging at T = T2. The re- 



sult eq.(3.2) combined with eq.(3.1) means that the 
mean domain size and the energy density are related 
by the scaling relation eq.( |2.6D . The latter implies 
that relaxation of the energy density in this regime 
is due to the coarsening of domain walls. 
The existence of isothermal regime 3), to which the 
system crossovers from regime 2) at around ^2 ~ ^mi''' i 
another important aspect of the aging dynamics observed 
by the present simulation. It is certainly a consequence 
of the persistence of domains through the T-shift process 
described before, and can be regarded as one of the mem- 
ory effects. An interesting question is if the phenomena 



in transient regime 2) are related with the rejuvenation, 
or chaos effect observed experimentally. 

3.4 Spm auto- correlation function 

We define the spin auto-correlation function in the T- 
shift process as 



[T2,Ti] 



t; tvji, t 



w2 7 ''wl J 



{Si{T + tml + ^w2)5'j(^wl + iw2))- 

(3.3) 

and measure relaxation of 

X[T2,T,]{T;t;t^i) = ^[1 - C[T2,Ti](-r;t2,^wi)] (3.4) 

with increasing ^2 after the T-shift. At each measure- 
ment, a fixed value of the time separation r is chosen as 
a parameter. 

We naively expect that the FDT holds for time regime 
r < tw2 as in the case of isothermal aging (see eq. ( 2 . 11 ) ) . 
If it is the case the above X[T2,Ti]{T]'t]'tvii) is identical to 
the ZFC susceptibility. In experiments, relaxation of the 
out-of-phase component of the ac susceptibility with a 
fixed frequency, say uj, is most frequently measured also 
in the T-shift protocol. It has been numerically ascer- 
tained in II that the relaxation of the ac susceptibility 
at a fixed frequency uj is very similar to that of the ZFC 
susceptibility with a fixed r = 27r/w = r^^. Thus we re- 
gard X[T2,Ti]{TLo','t','tvii) as essentially equivalent to the ac 
susceptibility measured in experiments. 

The function with tw2 = 0, denoteiLby Qtj.Ti] (''"; 0, twi) 
was examined in our previous worktS' hereafter referred 
to as III. We found that at r ^ follows sim- 

ilar scaling forms to those for CTiT'jt^) in isothermal 
aging described in §2.1 if the factor Lt (t) / Rt (tvi) in 
eq.(2T0) is replaced by i[T2,Ti](''";^wi)/-RTi(^wi)- Here 



L 



[T2,Ti] 



is a mean size of the droplet excitations. 



or quasi-domains, which are in local equilibrium at the 
new temperature T2 within time scale r after the T-shift. 

In Fig. 5 we show typical results of X[T2,Ti]{TLo','t','tvii) 
simulated for a fixed r^^ (= 32). Their behaviours are 
rather similar to those of ej^^ ^j] (t; twi) shown in Fig. 4. 
Just after the T-shift, it rapidly decreases (increases) 
and then relaxes slowly to the isothermal curve XT2 (tlo i) 
from below (above) for AT < (> 0). Furthermore, it 
also satisfies a relation similar to eq.(3.2): when the part 



(eff) 



it 



of its curve at ^2 > is plotted against t = t2 + t 
now relaxes to the isothermal curve from above (below) , 

X[T2,Ti]iTiv;t;twl) ~XT2(7"a. 



at t2^t^:f\ 



Here, a value of t 



(eff) 



(3.5) 

824) has been chosen so as to 



obtain a ^ood collapse of the two curves at the time range 



The resultant t,,,-, is again consistent with 



eq.( 3.12[) combined with eqs.(1.2) and (1.3) obtained in 



I. 



,(eff) 



The above feature in isothermal regime 3), ^2 ^ ^wi 
means that there the mean domain size and the cor- 
relation function are related by the scaling relations 
of eqs.(2.8), (2.9) and (2.10). Thus relaxation of 
X[T2,Ti]{TLo','t','tvii) in this regime is due to increase of the 
effective stiffness due to the coarsening of domain walls. 



Behaviour of 



X[T2,Ti](Tiv;t;twi) in transient regime 2) 
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Fig. 5. X[T2,Ti] ("^t^ 5 ^5 ^wi ) in a negative T-shift with Ti = 
0.8, T2 = 0.6 (a) and in a positive T-shift with Ti = 0.6, T2 = 
0.8 (b). The symbols and curves, including those in the in- 
sets, represent X[T2,Ti] ('^'^j ^wl) and XT{'^uj','t) in replacement 
of e|-7^2 ,Ti] (^5 ^wl ) and eqp{t) in Fig. 4), respectively. The values 
of twl and t^f^ are denoted in each figure. 



is proportional to AT = T2 — Ti . The boundary of the 
latter droplets may act as frozen-in domain walls and 
reduces the effective stiffness of droplets smaller than 
-^[T2,Ti](^2,^wi) just as in the case of the domain walls 



which separates different pure statesOO This yields 
the extra contribution, which is proportional to AT, to 
the relaxation. The scaling analysis on the excessive part 
in terms of _L[T2,Ti](^2, ^wi) and -R[T2,Ti](^; ^wi) is now un- 
der investigation. 

3.5 T -cycling process 

A T-cycling process is a combination of two T-shift 
processes with the same |AT| but with opposite signs. 
At a time ^2 = ^w2 after the T-shift from Ti to T2 the 
temperature is turned back to Ti , and the subsequent 
relaxation is observed. When tw2 ^ ^wf"'' T-change 
back to Ti is from an isothermally aged state at T2, 
as shown in Fig. 6(a). In this negative T-cycling pro- 
cess (T2 < Ti), the approximated ac susceptibility at 
^ = ^3 + ^w2 + ^wi, now denoted simply by XT{t){Tuj;t), 
merges into the isothermal curve xTi {tlo i) from below 



ts + ^l"*^'. Here 

t^'vf^^ is the effective waiting time given by ^w^^''''^'""^^'' 



when the former is plotted against t 
i"ecti' 

,(eff) 



(since tw2 ^ ^wi'^ ^^^^ case). In the case that tw2 
is small enough, we obtain a similar behaviour but with 
t^'vf^^ ~ twi (not shown). This means that aging at T2 
does not contribute at all to aging at Ti . This is nothing 
but the merapitjjfeffect and has been observed by many 
experiments JU'LStS' 

In a positive T-cycling process with T2 > Ti shown 
in Fig. 6(b), our simulated data still exhibit the mem- 
ory effect in the sense that the curve at ^3 > lies on 
the isothermal curve of Ti if it is shifted by a proper 
amount, i.e., t™*'''' which depends not only tw2 but also 
on twi- This type of memory retained over a positive 
T-cycling has been ipdeed observed in the ferromagnetic 
fine particle system.ES' 



is considered to refiect nature associated with quasi- 
domains. Let us here focus on the negative T-shift shown 
in Fig. 5(a). As pointed out just above, X[T2,Ti] ("'w ! ^5 ^wi) 
is definitely smaller than XT2 (''w at t ^ t^i . This is 
different from the rejuvenation, or chaos effect observed 
in experiments if it means that the system after the T- 
shift is younger than t^i as for an age in isothermal ag- 
ing at T = T2, i.e., X[T2,Ti] ("^o. ; ^; ^wi) as a function off 
with t'^fl— .smaller than twi coincides with XT2 (tlo i) at 
t2^trev.E3) It may suggest that overlap length /at intro- 
duced by the droplet theory and expressed as eq.(2.13) 
in §2.3, if it exists in the present SG model, should be 
larger than -Rti(^wi) observed in the present simulation. 
Let us then compare X[T2,Ti] ("'w ; i', ^wi) as a function of 

t'z with XT2 (''w ; ^"^wf'')- 'I'ii^ curves are shown in the 
inset of Fig. 5(a). Interestingly, the former is larger than 
the latter and the excessive part decays very slowly. This 
feature in transient regime 2) may be explained within 
the scenario sketched in §2.3 as the following. After 
^w2 from the T-shift, there remains excessive population 
of droplets of length scales larger than the size of the 
quasi-domain _L[T2_Tj](t2,^wi)- The excessive population 



3.6 Continuous T-change with an intermittent stop 

We have also performed simulations on a continuous 
T-change with an iutermittent stop. Following the ex- 
perimental protocol,E3) we cool a system by a constant 
rate from above Tc to below Tc with an intermittent stop 
at a certain temperature Tj below Tc, restart the cool- 
ing, and then heat it back by the same rate as in the 
cooling. A typical result is shown in Fig. 7. In this 
analysis XT{t){'''Lo',i) has been directly measured as the 
response to an ac field h[t) = hosm{uit) with hg = 0.3 
and 27r/w = 160 (MCS). 

At the intermittent stop at Tj, XT[t){'^'j^) relaxes 
downward and deviates from the reference curve of cool- 
ing. After restarting the cooling, XT(t)i'^''^) seen to 
merge into the reference curve. It appears that the long 
time aging at Tj does not affect subsequent cooling at 
temperatures lower than Tj by a certain amount. Here 
we call this phenomenon a rejuvenation-like effect, since 
Xy|,jj(w;t) merges with the reference curve but without 
a strong vertical increaa£.in contrast to the rejuvenation 
seen in the experiment .E3' 

When the system is heated back from the lowest tem- 





0.3 




0.28 




0.26 












0.24 


£ 








0.22 








0.2 




0.18 




0.16 



(a) 


\ 


?wi=512, rw2= 16384 
4=f«=2411 




71=0.8 


r,=o.6\ 




To,=27!/ft)=32 





Tatsuo KOMORI, Hajime YoSHINO and Hajime Takayama 

■2 



100 

(b) 



1000 ~ 

t,t 



10000 



100000 





0.3 




0.28 




0.26 












0.24 










0.22 








0.2 




0.18 




0.16 


Fig. 


6. X 


0.8 


, T2 


0.6 


, T2 



?wi=512, rw2=512 
4f«=5792 




71=0.6 



T.,=27i/ft)=32 



100 



1000 ~ 

t,t 



10000 



100000 



(efl) 



in Fig. 5 but with t = ts -\- 1 



perature of cooling process, XTmi'^j^) exhibits a dip at 
around T = T\, though it is rather shallow. This means 
that the memory imprinted by the long time aging at 
Tl in the cooling process, in fact, has not been erased 
until the system comes back to Ti from the lower tem- 
peratures. To conclude the present SG model exhibits 
the rejuvenation-like and memory effects within a time- 
window of the present simulation. 

§4. Discussions 

We have studied aging phenomena in various T- 
change protocols on the 3D Gaussian Ising SG model 
by Monte Carlo simulations. In particular, we have ar- 
gued that, in a T-shift protocol, slow aging dynamics 
after the T-change from Ti to T2 at ^2 = consists 
of two regimes: transient (^2 ^ ^wf'') ^^'^ isothermal 
(^2 ^ ^wf'') regimes, wh ere t i?^'' is the effective wait- 
ing time specified by eq.( ^.12j) . In the transient regime, 
re-distribution of thermal weights, from those at Ti to 
T2, of droplet excitations with the characteristic length 
scale i[T2,Ti](^2, ^wi) is taking place within each domain 
of the mean size i?[T2,Ti](^; ^wi)- We called the length 




Fig. 7. XT{t){^o:'^i) in a continuous T-change with an intermit- 
tent stop. We cool and re-heat the system continuously with 
a constant rate of J/3200 per MCS from T = 3.0 to T = 0.1. 
The open circles represent the reference curves in cooling (upper 
curve) and re-heating (lower curve) protocols. The solid circles 
represent the same protocol except for an intermittent stop in 
cooling at Ti = 0.8 for 32000 MCS. 



scale i[T2,Ti](^2, ^wi) as the size of quasi-domains. 



Suppose, in a negative T-shift process, we further de- 
crease temperature to T3 (< T2) at ^2 = ^w2 ^wf'') 
and wait a certain period. Then there appear now quasi- 
domains locally equilibrated to T3 within each quasi- 
domain of T2. Because of the wide separation of time 
scales at different length scales emphasized in §2.1, a 
nest (or hierarchy) of quasi-domains can be realized dy- 
namically when the system is cooled step-wise, or con- 
tinuously. At low temperatures the ac field can only 
excite droplets of the size corresponding to the 

frequency uj. If the size is small enough, a new contri- 
bution to the ac response will appear every time a new 
frozen-in domain wall appears associated with the quasi- 
domain due to a step of cooling. This is our interpre- 
tation of the rejuvenation-like effect seen in Fig. 7. Its 
observation then crucially depends on time and length 
scales by which one looks at the system. It is the easier 
to be observed, the shorter and smaller are the time and 
length scales. 

By means of our scenario the memory effect is rather 
easy to be understood. For example, the continuous T- 
change protocol shown in Fig. 7 is interpreted as follows. 
When the system is re-heated from the lowest tempera- 
ture, the hierarchical structure of quasi-domains is erased 
from the lower levels. When the temperature comes back 
to Ti, the large domains imprinted at the intermittent 
stop in the cooling process reappear and govern the re- 
sponse of the system, which causes the memory effect 
in this protocol. According to our scenario, a memory 
imprinted by aging at a certain temperature Tj is erased 
when the system is aged at different temperature, say 
Tj , for such a long period that the mean size of quasi- 
domains of Td catches up the domains imprinted at Tj , 
as what happens in regime 3) of the T-shift protocol. 
Even in that case, the memory of the previous thermal 
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history is kept in the effective waiting time, i.e., the sys- 
tem iooks more aged than the period that it is actuaiiy 
aged at T = Td (see eq.( ^.5| )). 

Let us now compare our simuiated resuits with exper- 
imentai observations. As pointed out aiready in §f , in 
the ferromagnetic fine particles system (FFPS), whose 
aging dynamics has been recently studied by Mamiya et 
al.f3 the microscopic time, Tmic, required for a mag- 
netic mom£iit of each fine particle to fiip is of the order 
of 10~^secE£r Therefore the time window of their obser- 
vation, 1 ~ 10^ ksec, in unit of Tmic is not much different 
from that of our MC simulation. It is also noted that 
in the FFPS each magnetic moment behaves as an Ising 
spin due to strong magnetic anisotropy. 

Indeed in the T-shift protocol, Mamiya et al. have 
observed relaxation of Xj^^ t; twi) which is quali- 

tatively quite similar to our simulated X[T2,Ti] (''o) ! ^5 ^wi) 
shown in Fig. 5. They have determined amount of the 
shift of X[y^ or t[^i \ by the same method 

that we have used for analysis of X[T2,Ti](''a) ! ^5 ^wi) in 
Fig. 5. The effective waiting time they obtained is 
written as t 
trast, t 



(eff) 



(eff) 



3 X (twi) with h ~ 1.0. By con- 

extracted from our simulation is given by 

eq.( 2.12| ), which, with making use of eq.(1.2), is reduced 
to tj^cx {t^if with h = z{T2)/z{Ti) (~ 1.11 (1.25) for 
Ti = 0.8 and T2 = 0.7 (0.6)). Also aging dynamics of the 
FFPS in the T-cycling protocol agrees qualitatively with 
the results of our simulation. In particular, the effective 
waiting time ti^^"^ which depends on t^i (memory effect) 
is observed even in the positive T-cycling process, though 
Mamiya et al have claimed that the memory is partially 
lost because some particles with a larger size than the 
average are affected by the chaos effect associated with 
the temperature change. 



For ordinary spin glasses with r^r 



10" 



-13g 

(eff) 



there 



have been only a few experiments, in which tj^j^ ' in the ill- 
shift process is systematically studied. Djurberg et alx3 
measured the ZFC magnetization of Cu:Mn 13.5at% spin 
glass with a fixed twi and varying AT. Via the FDT 
similar to eq.(2.11) their analysis is directly compara- 
ble with our simulation on C[T2,Ti](''"; 0, twi) in HI. The 
peak position of its relaxation rate with respect to Int 
may be assigned as t[^f^ (crossover from regime 2) to 
regime 3)). The experimentally observed t[^f^ becomes 
larger (smaller) than twi in the positive (negative) T- 
shift, which is qualitatively consistent with eq.( ^.12 ) with 
monotonic increase of l/z(T) on T. Quantitatively, how- 
ever, the shift observed experimentally is much larger 
than what is expected from eq.(1.3). 

Djurberg et alES> furthermore claimed that they ob- 
served the chaos effect in T-shift processes with |AT| 
larger than a certain value. The rejuvenation, or chaos 
effect in the T-shift and T-cycling protocols obs^ 
in Xy|,jj(w;t) has been reported for CdCri 7lno,3S4 
and Cu:Mn 2at%B spin glasses, as already^jpentioned 
in §1. The observed XT(t)(^''^) claimedO-' to indi- 
cate that the system is literally rejuvenated in the sense 
described in §3.4. Neither such a chaos effect nor lit- 
eral rejuvenation have been_seen in our present simula- 
tions (see also Kisker et alx3). At the moment it is not 



ml 



clear yet whether this discrepancy originates simply from 
the difference in time scales of the observations, or from 
something else. 

To conclude, we have simulated various T-change pro- 
tocols of aging in the 3D Ising SG model, and have pro- 
posed a scenario in terms of domains and quasi-domains: 
the former ever grow in any aging process in the SG phase 
and the latter appear within the former or within them- 
selves. The results observed in our simulations are rea- 
sonably well interpreted by the scenario. Qualitatively, 
it is also consistent with many features observed exper- 
imentally, particularly, those related to the memory ef- 
fect. For certain phenomena such as the growth law of 
Rxit), the consistency is even semi-quantitative. Con- 
cerned with the rejuvenation, or chaos effect, however, 
the scenario may not be applicable, and further studies 
are required. 
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